Measurement of the infrared transmission through a single doped GaAs quantum well 
in an external magnetic field: Evidence for polaron effects 
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Precise absolute far-infra-red magneto-transmission experiments have been performed in magnetic 
fields up to 33 T on a series of single GaAs quantum wells doped at different levels. The transmission 
spectra have been simulated with a multilayer dielectric model. The imaginary part of the optical 
response function which reveals new singular features related to the electron-phonon interactions has 
been extracted. In addition to the expected polaronic effects due to the longitudinal optical (LO) 
phonon of GaAs, a new kind of carrier concentration dependent interaction with interface phonons 
is observed. A simple physical model is used to try to quantify these interactions and explore their 
origin. 



Polaronic effects, due the Frohlich interaction between 
electrons and the longitudinal optical (LO) phonon of a 
polar semiconductor, have been the object of many re- 
ports [3, 0, In quasi- two-dimensional (Q2D) GaAs 
based structures, there is evidence of free polaronic ef- 
fects (PE) but few experiments 0, [B| have focused on 
cyclotron resonance (CR) in the region above the Rest- 
strahlen Band energy (RBE) of GaAs (this requires for 
GaAs a magnetic field strength beyond 23T). Theoreti- 
cally PE were first studied by Lee and Pines [6(] and later 
by Feynman [7| . It was realized that this kind of interac- 
tion could not be properly handled by perturbation the- 
ory and requires a global treatment. Such an approach 
was proposed by Feynman et al. 8], and referred to as 
the FHIP model. It has been invoked to explain the po- 
laronic mass and later extended by Peeters and Devreese 
Q to extract the conductivity of the Q2D electron gas in 
the presence of PE. In the latter case, its major effects 
are expected to be observed in the imaginary part of the 
response function at energies larger than the RBE. This 
prediction has motivated the present study. 

Far infra-red magneto-optical experiments have been 
performed at magnetic field strengths up to 33 T and 
at a fixed temperature of 1.8 K, on a scries of single 
modulation-doped quantum wells (QW) of width L = 
13 nm with different doping levels N s ranging from 2 to 
7.7 x 10 11 cm -2 and mobility exceeding 10 6 cm 2 /(V.s). 
The structure of these samples is similar to those re- 
ported earlier a single GaAs QW is sandwiched be- 
tween two GaAs-AlAs superlattices 5-doped symmetri- 
cally with Si-n type dopant on both sides of the QW. 
In the present case however the epilayer is not lifted-off 
from the GaAs substrate. As a result the samples are 
optically opaque in the RBE range. For each fixed value 
of the magnetic field B, an absolute magneto infra-red 
transmission spectrum (TA(B,u>)) is measured by using 
a rotating sample holder containing a hole to obtain a 
reference spectrum under the same conditions as for the 



sample. The Faraday configuration is used with the k 
vector of the incoming light parallel to B and also to the 
growth axis (100) of the sample. These spectra TA(B, u) 
are in turn divided by TA(0, ui) to obtain the relative 
transmission spectra TR(B,lj) which will be displayed 
in the present paper. The analysis of these spectra is 
based on a multi-layer dielectric model [l(|. This is es- 
sential because, in the frequency range of interest, the 
spectra can be distorted by dielectric interference effects 
even in the absence of electron-phonon interactions. 

The TR spectra of two characteristic samples, sample 
MA-2226 (N s = 3 x 10 11 cm" 2 ) and sample MA-1490 
(N s = 6.3 x 10 11 cm -2 ), are displayed in Figs.QJ, andQJl 
respectively. For sample MA-2226, the CR absorption al- 
most vanishes just above the RBE for B ~ 24 T. But as 
B is increased beyond 24 T it starts to increase and also 
its lineshape changes with field. A " singularity" appears 
to occur when the CR energy approaches 45 meV (the 
transverse optical phonon (TO) energy of AlAs as mea- 
sured in the TA(0) spectrum) at B ~ 27 T. Its linewidth 
peaks when its energy is around 48 meV (B ~ 29 T) be- 
fore recovering its original line shape below the RBE at 
even higher fields. In contrast, for sample MA-1490, the 
CR transition is clearly seen when it first emerges above 
the RBE, its lineshape has broadened already when its 
energy is around 38 meV (B ~ 24.5 T), it shows some sin- 
gular behavior at 45 meV (B ~ 28 T) while its linewidth 
appears to reach a maximum at the same time. Finally 
it recovers the original low fields lineshape for B > 32 T. 

To simulate these complex spectra, we start with the 
dielectric function e of the doped QW in addition to the 
appropriate dielectric functions of undoped barrier layers 
[lfj . For the doped QW, the diagonal part e xx is written 
as: 



uj[uj - (lo np - Re(£)) + i(r) + Im(S)] 
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FIG. 1: (Color on line) (a) and 
(d): Relative transmission spectra 
of samples MA-2226 and MA-1490 
respectively, for different magnetic 
fields, (b) and (e): simulation of 
some experimental spectra (open 
dots) with the multidielectric model 
( continuous curve): each curve is 
shifted by 0.5 with a scale enlarged 
as noted, (c) and (f ) Fitted parame- 
ters hu) c /B and So as extracted from 
Eq. (1) (open dots) for both sam- 
ples; the dashed lines in the upper 
panels (So) are the resulting decom- 
position of the polaronic and XP in- 
teractions adopted to fit Im(£); the 
total contributions of the interac- 
tion (continuous curves in (c) and 
(f)) are used to simulate the spec- 
tra in (b) and (e) respectively; the 
oblique dashed lines in the lower 
panels of (c) and (f) mimic the 
non-parabolicity contribution with 
the corresponding energy scale (see 
text) 



in which the x-axis is assumed to lie in the plane of the 
QW, £l is the contribution of the GaAs lattice, lonp 
is the contribution to the CR frequency of the non- 
parabolicity (NP) effects [llj, rj is the field independent 
contribution of the background defects to the line width 
while represents the self energy correction due to 

interactions between the electron gas and elementary ex- 
citations to be discussed further. The effective mass m* 
entering the plasma frequency io p is calculated with the 
same NP model. We note that the lineshape and inten- 
sity of the CR transitions at low and high fields are prac- 
tically the same. Hence we can assume that there is no 
loss of carrier density during the magnetic field runs. In 
the fitting process, we first neglect the frequency depen- 
dence of S and are left with two independent parameters 
lu c = ujnp — Re(£) and Sq = rj + Im(E) for each mag- 
netic field. The resultant parameters (hio c /B and Sq) are 
displayed in Figs. [TJ; and QJ (open dots) as a function 
of the energy Tilo c (B) (in these figures some typical er- 
ror bars reflect the estimated uncertainties of the fitting 
procedure). From these plots we notice that below the 
RBE, <5n is constant and weak being determined mainly 
by rj. Above the RBE there are sudden and large in- 
creases in 5q for both samples. However their behaviors 
are quite different. For sample MA-2226, there are two 
distinct contributions to Sq: one lower field component 
starts very large but decreases strongly with field and 
a second component which starts from almost zero but 
goes through a maximum at B corresponding to about 
48 meV. For the time being we will label the interaction 
which contributes to this peak in <5n as XP. For sample 



MA-1490, only the XP contribution to Sq is observed. In 
this case the peak occurs at B corresponding to 42 meV. 

Thus our results suggest that there are two kinds of ex- 
citations interacting with the Q2D electron gas at high 
fields: one which is observed in the lower doped sam- 
ples and another one labelled as XP which is observed 
in both samples. Very similar results have been ob- 
tained with another sample MA-2227 doped at the level 
N s — 2 x 10 11 cm~ 2 . In this sample the first interaction 
is a little bit stronger than for sample MA-2226 while 
the XP peak is quite similar. We assign the first interac- 
tion to polaronic effects since it appears to peak at fields 
when the CR frequency will resonate with the GaAs LO 
phonon frequency (slab mode). This assignment is con- 
sistent with the observation that this interaction is very 
sensitive to the carrier concentration. In sample MA- 
1490 this interaction seems to have disappeared from So- 
We note that in a previous report [H[ , on lift-off samples 
more highly doped in the range N s = 7 to 9 x 10 11 cm~ 2 , 
only some weak interaction around the TO of GaAs was 
observed. This result was later on analyzed by Klimin 
et al. [l§] as PE with the LO phonon screened by elec- 
trons to the point that its value becomes close to the TO 
phonon frequency. The same analysis can be applied to 
the sample MA-1490 and other higher doped samples. 

To mimic PE we restrict ourselves to one polaron first 
and, using the FHIP model, write down Im(S(w)) as 

Im(£H) = u LO —aF\A Q {u)\ l l 2 e- R \ A ^Q{Mu)) 

LO 

(2) 

where Aq(lj) = lj/loq — 1, 0(x) = 1 for x > 1 or zero 
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otherwise. In the absence of screening ujq — ujlo- If 
ujq is known, Eq. (2) depends on two parameters aF 
and R. The FHIP model, developed in a one electron 
picture, assumes that electron-phonon interaction is har- 
monic instead of Coulombic as in the case of the Frohlich 
interaction. The model depends on two parameters: v 
and w in reduced units of ujlo- One of them (w) is 
close to 1 while v 2 — v? represents the force constant of 
the harmonic interaction (equivalent to the Frohlich con- 
stant a). In this case R = (i? / v? — l)/v whereas aF is 
a global quantity which depends on v and w but should 
also imply corrections for the dimensionality of the prob- 
lem and for screening effects. The present approach can 
be regarded as a test of the FHIP model when applied to 
Q2D electrons under high magnetic fields. 

Using this simplified model, one can fit the data on 5q 
for the polaronic contribution. For sample MA-2226 we 
have to add contribution for the XP interaction (Fig. [TJ;, 
upper panel) . If we assume that the XP interaction sat- 
isfies the linear response function theory, then Im(£(cj)) 
and Re(£(u)) are related by the Kramers-Kroning (KK) 
relations. The resultant constraint on the fitting process 
is that the KK transformation of Sq(u>) should reproduce 
the variation of TujJ c (lo)/B. One has also to fit the NP ef- 
fects, h(VNp(B)/B versus ?iujnp(B) with standard mod- 
els [llj (dashed oblique lines in Fig. [TJ:) but this corre- 
sponds simply to a global shift of the curves. The fitted 
functions hui c (ui)/B and Sq(uj) are then inserted in the 
multilayer dielectric model from which one can compute, 
for each value of B, the TR spectrum and compare it to 
the corresponding experimental spectrum (Fig. [TJd) . The 
agreement is quite satisfactory as one can see from that 
figure. However, it is clear that the fitting of PE is not 
unique for these reasons: (i) we have assumed for sample 
MA-2226 that ujq = ^lo neglecting, therefore, screening 
effects, (ii) we have treated the problem with one polaron 
instead of several polarons. Even if this effect does not 
influence Im(S(cj)), the KK transformation depends on 
all contributions to Im(E(w)) including those at higher 
energies than we can reach experimentally. The same 
procedure has been applied to simulate the data on sam- 
ple MA-1490, but here we have assumed a weakened PE 
starting near the TO phonon energy of GaAs in order to 
reproduce the non-linearity of T%lo c (lo)/ B below the RBE 
(lower panel of Fig. QJ) . The choice of the corresponding 
fitting parameters in Eq. (2) is now completely arbitrary. 
We can also reproduce quite well the experimental trans- 
mission spectra as shown in Fig. [T]3. 

Although we do not know the value of the pre-factor 
aF in Eq. (2), it is instructive to compare the parame- 
ter R entering this equation. One finds R « 25 and 18 
for samples MA-2227 and MA-2226 respectively. These 
values could, of course, be lowered if one assumes some 
screening of the LO phonon in Eq. (2) and if we in- 
clude higher polarons in the analysis but they remain, at 
least, an order of magnitude higher than that proposed 
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FIG. 2: (Color on line) Variation of Im(E) with the filling fac- 
tor v for different samples: MA-2227 with open dots (7V S = 
2 x 10 11 cm" 2 ), MA-2226 with full dots (N B = 3 x 10 11 cm" 2 ), 
MA-1490 with open squares (N s = 6.3x 10 11 cm" 2 ), MA-1523 
with full triangles (N s = 6.6 x 10 11 cm" 2 ), MA-1315 with 
open stars (N a = 7.1 x 10 11 cm -2 ), MA-1316 with full squares 
(N s — 7.7 x 10 11 cm -2 ). Insert: Variation of the character- 
istic energies of the maximum of the XP interaction 
with the corresponding magnetic field i? ma a;(see text). The 
corresponding error bars take into account the non symmetric 
shape of the interaction. 



(i? = 0.04) for unscreened PE [13j. These values imply 
that the v and w parameters in the FHIP model have 
to be smaller than 1 even for the samples MA-2226 and 
MA-2227. This raises the question of the physical inter- 
pretation of these parameters: knowing that ujlo cannot 
be smaller than uito, it seems that a more physical ap- 
proach should be to relate in some way these parameters 
to the TO-LO splitting rather than to a single frequency 
like ujlo- 

We focus next on the discussion on the possible origins 
of the XP interaction. The variation of Im(S) = <5o — rj 
is plotted in Fig. [5] for six different samples as a function 
of the filling factor v = N s Qq/B where $o is the flux 
quantum. The XP interaction which is only observed for 
energies larger than the RBE, clearly decreases for the 
higher doped samples (MA-1315 and MA-1316) and in- 
deed disappears completely when N s > 9 x 10 11 cm -2 . 
This result indicates that the XP interaction is strongly 
screened by free electrons. The energy E^ x correspond- 
ing to the maximum of this interaction decreases while 
the corresponding value v max , the value of v where this 
maximum occurs, increases or E^ x increases with the 
corresponding value of B max (insert of Fig. [2]) . This inter- 
action is clearly dissipative and in these samples of very 
high mobility, is unlikely to be caused by impurities. We 
are then left with intrinsic mechanisms such as electron- 
phonon interaction as a plausible explanation knowing 
that in these samples the inter-subband energy is of the 
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order of 65 meV and cannot play a role. In a GaAs 

QW, sandwiched between AlAs layers, several kinds of 
optical phonons are present. The most well known ones 
are the confined phonons (slab modes) mainly associ- 
ated with mechanical motion of atoms. The frequency 
of these modes lies inside the RBE. The other phonons 
of dielectric origin (which could be hence screened by 
free carriers) are the interface phonons. Because of the 
reflection symmetry of the GaAs QW with respect to its 
center, these modes are divided into either symmetric or 
anti-symmetric modes [l4T | . The symmetric modes which 
are infra-red active, have been invoked in Ref. Q to ex- 
plain the "splitting" of the CR transition observed near 
fouiTO {AlAs) (45 meV). Such splitting is observed also in 
all our samples including the very heavily doped sam- 
ples (up to N s — 1.9 x 10 12 cm -2 ). As seen in Figs, [lb 
andQJ, it can be perfectly reproduced by our simulation 
with the multi-layer dielectric model as a result of pure 
interference effects in the transmission spectra. We also 
note that the wave function of these symmetric interface 
modes does not vanish at the center of the QW (unlike 
the anti-symmetric modes) and therefore have significant 
overlap with the wave function of the Q2D gas. Thus the 
idea of possible interaction between the interface modes 
and the Q2D electron gas deserves more detailed analy- 
sis. The solutions for these modes between a GaAs layer 
sandwiched between two AlAs layers are given by [P4j : 



definitive understanding into all the interactions between 
the Q2D electrons and other elementary excitations like 
phonons, one should obtain experimental results on lift- 
off samples which can provide more unambiguous values 
of the threshold energies related to the simple polaronic 
effect. 

In conclusion, we have performed infra-red magneto- 
optical transmission measurements of a Q2D electron gas 
in a single modulation-doped GaAs QW with different 
densities, for magnetic fields high enough to scan the 
cyclotron resonance frequency beyond the Reststrahlcn 
band of GaAs. From the experimental spectra, we have 
extracted the imaginary part of the response function 
which reveals several singularities whose number and 
strength depend on the carrier density. These singulari- 
ties have been attributed to electron-phonon interactions. 
One of these interactions involving the LO phonon of 
GaAs has been treated with a simplified version of the 
FHIP polaronic model. But to explain all our results 
quantitatively, it is necessary to invoke a more elaborate 
theory which includes effects of screening and possible 
interaction with interface modes. 
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SAlAs = SGaAs x tanh(q//L/2) (3) 

where qi / is the wave vector of the interface electromag- 
netic wave travelling parallel to the plane of the Q2D 
gas. In the absence of magnetic field, the two solutions 
of Eq. (3) lie in the Reststrahlen bands of GaAs and AlAs. 
But, when the magnetic field effect is included in EcaAs, 
through Eq. (1) for instance, the results are different: (i) 
depending on the parameters the number of solutions of 
Eq. (3) can be larger than 2, (ii) some of the solutions 
are now strongly dependent on the magnetic field and 
can extend to energies higher than Hluto(AIAs) as ob- 
served here. When increasing N s , one would expect that 
an interaction with the Q2D electron gas will result in a 
downshift (or renormalization) of the interface mode en- 
ergy in a way similar to what is observed for the polaronic 
effects in Figs. [I}; andQJ. This effect has never been inves- 
tigated to our knowledge but certainly deserves further 
studies. In particular it is desirable to understand the 
mechanism of clcctron-phonon interaction which makes 
this interaction resonant for a specific energy E^ x . 

At present we prefer the explanation of the XP in- 
teraction in terms of interface modes though we cannot 
exclude more speculative interpretations, such as, strong 
non-linear effects of the polaronic interaction. 

It is clear that in order to get a deeper and more 
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